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Abstract 



In this paper, we use various anstazes motivated from our earlier works on transonic gas flows, 
boundary layer problems and Navier-Stokes equations to find new explicit exact solutions with multiple 
parameter functions for the equation of geopotential forecast and the equations of nonlinear magnetohy- 
drodynamics. 

1 Introduction 

Partial differential equation 

{HxX + Hyy)t + Hx{Hxx + Hyy)y " H y{H xX + H yy) X = ^Hx (l-l) 

is used in earth sciences for geopotential forecast on a middle level (e.g., cf [Ki] and Page 222 in 
[I]), where A; is a real constant. Kibel' [Ki] found the Gaurvitz solution of the above equation. 
The well known Syono solution was given in [Sy]. Katkov [Kal, Ka2] determined the Lie point 
symmetries and obtained certain invariant solutions. The other known solutions are related to 
the physical backgrounds such as configuration of type of narrow gullies and crests, flows of type 
of isolate whirlwinds, stream flow, springs and drains, hyperbolic points, and cyclone formation 
(e.g., cf. Pages 225, 226 in [I]). 

In magnetohydrodynamics, it is very important to study the nonlinear MHD equations: 

Ipt + Vxtpy - ^ytpx = -^z, (1-2) 



{y^xx + ^yy)t + ^xiy^xx + Vyy)y " Vyi'fxx + ^yy)x 
= {'^xx+1pyy)z + 1px{'>Pxx+1pyy)y-^y{'4^xx+i^yy)x, (1-3) 

where i-p and ^ip are the potentials for the velocity and the transverse component of the magnetic 
field, respectively (they can also be interpreted as the potential of an electric field and the z- 
component of the vector potential of the magnetic field) (e.g., cf. [KP] and Page 390 in [I]). 
We refer [P] for more information on magnetohydrodynamics. Samokhin [S] (1985) determined 
the Lie point symmetries of the above equations, conservation laws and some solutions in terms 
the solutions of the other partial differential equations. Bershadskii [B] found a connection 
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between the energy conservation law and the uniqueness of classical solution of the nonlinear 
MHD equations. The reason of solving the geopotential equation together with the nonlinear 
MHD equations is that the left hand side of the equation (1.1) coincides with that of the equation 
(1.3). 

Based our earlier works on transonic gas flows [XI] , boundary layer problems [X2] and Navier- 
Stokes equations [X3], we give in this paper various ansatzes related to algebraic characteristics 
of the above equations to find new explicit exact solutions with multiple parameter functions. 
By specifying these parameter functions, one can obtain the solutions of certain initial-value 
problems of the above equations. 

The symmetry group of the geopotential equation (1.1) is generated by the following trans- 
formations: 

Ta,b{H) = H{t + a,x,y + b), Tc{H) = c^H{c-H,cx,cy), a,b,c e R, c ^ 0; (1.4) 



T^,p{H)=H{t,x + a,y) + a'y + P, (1.5) 

where a and /3 are arbitrary functions of t. The symmetry transformations of the nonlinear 
MHD equations (1.2) and (1.3) that we are concerned with are: 

Ta,b{'P) ='P{t + a,x,y,z + b), Ta,b{tp) = il^it + a,x,y, z + b), a,beM; (1.6) 

Ticiif) = c" V(ct, X, y, cz), Ti^(^) = V(ct, x, y, cz), (1.7) 

T2c{y^) = c~^ip{t, cx, xy, c), T2c{'ip) = c~^i^{t, cx, cy, z) (1.8) 

for ^ c G M; 

Ta,T{v) =^{t,x + a,y,z)+aty + Tt, Ta,r{''P) = '(p{t,x + a,y, z) + a^y + t^, (1.9) 

Taiv) = v{t,x,y + a,z) -atx, Ta{'(p) = '(p{t,x,y + a, z) - a^x, (1.10) 
where a and r are any functions of t,z; 

Taijp) = ^{t,xcos la -I- y sin 2a, —a; sin 2a -|- y cos2a, z^ -|- ol{x^ + y^), (l-H) 

Taifp) = ijjitjX cos 2a + y sin 2a, —x sin 2a -|- y cos 2a, z) + a'{x^ + y^), (1-12) 

Tp{ip) = ip{t, X cos 2/3 + y sin 2/3, -x sin 2P + y cos 2/3, z) + /3'(x^ -|- y^), (1.13) 

Tpitp) = 'tp{t, X cos 2/3 + y sin 2/3, -x sin 2p + y cos 2/3, z) - I3'{x^ + y^), (1.16) 



where a = a{t + z) and /? = f3{t — z) are arbitrary one-variable functions. The above transfor- 
mations change solutions to solutions. For convenience, we always assume that all the involved 
partial derivatives of related functions always exist and we can change orders of taking partial 
derivatives. We also use prime ' to denote the derivative of any one-variable function. 

In Section 2, we solve the equation (1.1) of geopotential forecast. We find explicit exact 
solutions of the nonlinear MHD equations (1.2) and (1.3) in Section 3. 

2 Solutions of Geopotential Forecast Equation 

In this section, we find two families of exact solutions of the geopotential forecast equation (1.1). 
Let a and (3 be functions of t. Set 

w = ax + Py. (2.1) 

Assume 

H = <i){t,w)+ /ly^ + TX + vy, (2.2) 



where ^ is a two- variable function and r, /n, v are functions in t. Note 

H^ = a4>^ + T, Hy = p(f)^ + 2iiy + u, H^^ + Hyy = 2/1 + {a^ + p'^)(f>^^, (2.3) 
{H^^ + Hyy)t = 2fi' + (a^ + + + + {a'x + (2.4) 

{H^^ + = + (^^o^cc + i^yy)?; = (a^ + P'^)P(l>^u,^- (2.5) 

Thus (1.1) becomes 

2/x' + (a^ + /-J^)',^^^ + (a^ + /32)04^^ - A;(a(/.^ + r) 

+(a2 + + (/?' - 2a//)y + /?r - ai/]<^,:^,^^ = 0. (2.6) 

In order to solve the above equation, we assume 

2ii' = kT, T = ad", u = l3d'', (2.7) 

for some function i? of t, and 

a'x + {0 - 2an)y = 0. (2.8) 
Note that (2.8) is equivalent to the following system of ordinary differential equations: 

a' = 0, P' - 2aii = 0. (2.9) 

By the first equation in (2.9), we have a = c G M. So r = cd'' according to the second 
equation in (2.7). Moreover, the first equation in (2.7) yield 

;. = ^^, coGM. (2.10) 
Hence the second equation in (2.9) becomes 

fi' - c{kcd' + cq) = {). (2.11) 

Therefore, 

13 = c{kc-d + Cot) + d, deR. (2.12) 
According to the third equation in (2.7), 

I' = [cikcd + Cot) + dj'd". (2.13) 

Now (2.6) becomes 

(a^ + P^cp^^ + (a^ + /?2)0i^^ - kce'^(l)^ = 0. (2.14) 
Modulo the transformation in (1.5), it is enough to solve the following equation: 

(a^ + py<p^ + (a^ + f3'^)4>t^ - kce^4> = 0. (2.15) 
The above equation can written as 

[{a^ + p'^)(f>^]t - kc(j) = 0. (2.16) 

So we take the form 



(l){t,-Uu) (f){t,Tu) 



a2 + /32 c2 + [c(fcci? + Co*) + df ' 



(2.17) 



Then (2.16) becomes 

" c2 + [c(/cC1? + Cot) + d]2" 

Thus we have the solution: 



(2.18) 



^ = g 7 0^ + [c(^c^ + cot) + d]2 + '^^^J 

X sm ^6.^ + - ^ y c^ + dV^c^^ J ' (2-19) 
where a^, ftj, Cj, dj are real constants such that (a^, bi) / (0, 0). 

Theorem 2.1. Let t9 6e any function oft and let ai,bi,Ci,di,Co,c,d for i = 1, ...,m be real 
constants such that {c,d),{ai,bi) 7^ (0,0). We have the following solution of the geopotential 
forecast equation (1-1): 

H = fec^' ^''% ^ + ^"[cx + [cikcd + cot) + d]y] + ^ 



2 ^ ' c2 + [c(A;ci? + cot)+d]2 

(fecft- f dt \ 

bi[cx + [c{kc-& + cot) + d]y] +Ci- — ^ / — ^ . (2.20) 

^ / jyj af + bjj c^ + [c{kc'd + Cot) + df J ^ ' 

Next we set 

w = x^^y^. (2.21) 

Assume 

H = ^{w) - y (2.22) 

where ^ is a one-variable function. Note 

H^ = 2x^', Hy = 2y^'-\, H^^ + Hyy = 4{^' + w^"), (2.23) 

{H^^ + Hyy)^ = 8x{2C" + w^"'), {H^^ + Hyy)y = 8y{2^" + wC"'). (2.24) 
Then (1.1) is equivalent to: 

A{2C" + we") = k^'. (2.25) 
Modulo the transformation (1.5), we only need to solve the equation: 

e + we' = ^^. (2.26) 

To solve the above ordinary differential equation, we assume 

00 

C = ^ro^aj + 6ilnro), ai,bieR. (2.27) 

i=0 

Observe 

00 

= '^w'-^{iai + bi + ibilnm), (2.28) 

i=0 

00 

= - + - ^)bi + i{i - l)6ilntu). (2.29) 



So (2.26) becomes 



oo 



k 

'^w'~^{i'^ai + 2ibi + i^bihiw) = - ^w\ai + bilnw), (2.30) 



equivalently, 
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i=0 i=Q 



k k 

{i + ifai+i + 2(i + l)hi+i = -Ui, {i + 1)^6^+1 = -hi. (2.31) 



Hence 

_ bpk' _ gpfc' \^ rr, or,N 

Thus 

Theorem 2.2. Zei 6 and c 6e any reaZ constants. We have the following steady solution of 
the geopotential forecast equation (1-1): 



Remark 2.3. Applying the transformation (1.5), we obtain the following non-steady solu- 
tion: 



j=0 

^2 , „.2\ 



(i!)24^ 

where a and (3 are arbitrary functions of t. 

3 Solutions of Nonlinear MHD Equations 

In this section, we will find multiple parameter function exact solutions of the nonlinear MHD 
equations (1.2) and (1.3). 
We first assume 

(p = atxy + fx + gy + h, tp = a^xy + Xx + fiy + p, (3.1) 

where f,g,h,a,X,iJ, and p are functions in t,z. Then (1.3) naturally holds. Moreover, (1.2) 
becomes 

(At - fz)x+ {pt - gz)y + Pt - hz + {(Tty + /) (cr^x + p)- {atx + g) {a^y + X) = 0, (3.2) 
equivalently, 

Xt-fz + (Tzf - atX = 0, (3.3) 

Pt- gz + cftp - cFzg = 0, (3.4) 

Pt - hz + fp- gX = 0. (3.5) 



Solving (3.3) and (3.4), we get 

f = ^te'', 9 = ne-'', A = i?,e'^, /x = r^e"'^, (3.6) 
where ■& and r are arbitrary functions in t, z. Moreover, (3.5) becomes 

Pt-h, + ^tTz - ^zTt = 0. (3.7) 

Thus 

h = -dtr, p = ^zT (3.8) 
modulo the transformation of type Tq^^ in (1.9). So we have the following simple conclusion. 

Proposition 3.1. Let cr, and r he functions oft,z. We get the following solution of the 
nonlinear MHD equations (1-2) and (1-3): 

^ = atxy + ^te'^x + Tte^'^y + ^tT, (3.9) 

tp = azXy + -&ze'^x + Tze''''y + -&zT. (3.10) 

The main objective in this section is to find more sophisticated exact solutions of the equa- 
tions (1.2) and (1.3). Let 9 and e be functions in t,z. Set 

w = Qx + ey. (3.11) 

Suppose 

ip = Cit,z,w) + atxy, ip = r]{t,z,w) + azxy, (3.12) 

where ( and r] are functions in t,z,w to be determined, and a is a function of t,z. Then (1.2) 
becomes 

Vt-Cz + + ety - ati^x - ey)]ri^ - [SJ^x + e^y - a^i^x - ey)]C-c^ = (3.13) 
and (1.3) becomes 

(9^ + e^)tC^^ - (9^ + e'^)zV^^ + (9^ + e^){Ctz,^ + [9tx + ety - ati^x - ey)]C^^^ 

-■nzvjw - [^zx + Ezy - az{^x - ey)]r]^^^} = 0. (3.14) 

In order to solve the above system of partial differential equation, we assume 

$JsX + £sy - (7s(9a; — ey) = 0, s = t,z, (3.15) 

equivalently, 

^5 — (TsQ = 0, ffs + CTsE = 0, s = t,z. (3.16) 

So 

cj = 6e'', e = ce-^, 6,cGM. (3.17) 

Moreover, (3.13) becomes 

r}t-Cz = 0. (3.18) 

Hence 

Q = Ft{t,z,w), ri = F,{t,z,w) (3.19) 



for some three variable function F. Now (3.14) becomes: 

+ e^Ft^^ - (9^ + e^),F,^^ + (9^ + e'){Ftt^^ - F,,^^) = 0. (3.20) 
Modulo the transformation in (1.9), we have 

F = F + \w (3.21) 
for a function A of t, z and a function F of t, z, w such that 

($>2 + e\Ft - (^2 + e'),F, + (9^ + e'){Ftt - F,,) = 0. (3.22) 

Rewrite (3.22) as 

In order to find exact solutions of the above equation, we take the following special cases of 
9^ + e'^. Without loss of generality, we assume 6 = 1 in (3.17). 

Case 1. 3=2 + £2 ^ ^a^t+a^z ^i^j^ oi, 02 G M. 

Our assumptions says 

e^" + c^e-^^ = e"l*+"2^ (3.24) 

(e2<^)2 _ gait+a22g2a + ^ = q. (3.25) 

e^'' = — (3.26) 

or 



equivalently. 



So 



^ gOit+agz _ ^g2(ait+a22) _ 4^2 



Hence 



or 



= . (3.27) 

a = i (|ln (|e"i*+"22 + yg2(ait+a2z) _ 4c2^ _ In 2) (3.28) 



cr = ^ (^In (^e"i*+"2z - Ve2(«i*+»2^) - 40^) - ln2) . (3.29) 
Now (3.23) is equivalent to 

Ftt + aiFt = F,, + a2F,. (3.30) 

For ^ a G M, we denote 

V{a, s) = ads + (3.31) 

and 

6,o(a, s) = 1, Ci,i(a, s) = Y: tSCl ^ (3-32) 

6,o(a,s) = e— , ^2,^(a,5) = (-l)»e""-|] (3.33) 
for < z G Z. Moreover, we let 

2i 2i+l 

V{0,s) = dl ^,.(o,s) = j^, ^2,i(o,s) = ^2Hri)! ^^■^^'^ 



for < i € Z. Then 



or 



V{a, s)(C.,o(a, s)) = 0, V{a, s){^e,i{(^, s)) = Ce,i-i{a, s), e = 1, 2, (3.35) 
for < z G Z. Given one-variable functions 

{ae,i,{w),l3e,j,{w) I e = 1,2; = 0,1,..., m-e; je = 0,1,..., Ue}, (3.36) 
we have the following solution of (3.30): 

me ne 

^ = XI [^(^^A^)^eA(^l^''^)^hme-iia2,z) +^Pe,j{^)^e,jiai,t)^2,ne-jia2,z)]. (3.37) 
e=l,2 i=0 j=0 

Case 2. Case 1. + e"^ = e°-^^z°-^ with ai, a2 G M such that 02 {0, Z + 1/2}. 
As (3.24)-(3.28), we have: 

(7 = ^ (in (^e"i*z"2 + Ve2«i*z2«2 - 4c2) - In 2) (3.38) 

a = ^ (^In - Ve2»i*z2«2 - 4c2) - In 2) . (3.39) 

In this case, (3.23) is equivalent to 

Ftt + aiFt = F^, + a2Z-'^F,. (3.40) 

For a G M \ {0,Z + 1/2}, 

Ci,o(a, s) = 1, Ci,i{a, s) = -, (3.41) 

^ '''■[[r=ii^ + 2r - 1) 

2i+l-a 

C2,o(a, .) = C2,^{a, s) = (3.42) 

2*z!ll^^i(2r + 1 - a) 

for < i G Z. Denote 

P(a,s) = 92 + ^9,. (3.43) 

Then 

Via, s){Ce,o{a, s)) = 0, Via, s){Q,i{a, s)) = Ce,i-i{a, s), e = 1, 2, (3.44) 

for < i G Z. Given the one- variable functions in (3.36), we have the following solution of 
(3.40): 

me rie 
^ = XI i^^'^A^)^e,i{ai^iKi,me-i{a2,z) + ^(5e,j{^)^e,j{ai,t)C2,ne-j{a2,z)]. (3.45) 
e=l,2 i=0 j=0 

Case 3. Case 1. 'ks'^ + e"^ = t^iz"^ ^j^h ai, 02 G M \ {0, Z 1/2}. 
As (3.24)-(3.28), we have: 

cr = ^ {\n (t"'^z"-'' + \/i^«iz2a2 _4c2^ _ In 2) (3.46) 



or 

(7 



^ (in (^e^z"^ - ^t2ai^2a2 _4c2^ - In 2^ . (3.47) 



In this case, (3.23) is equivalent to 

Fu + a^t-^Ft = F,, + a2Z-^F,. (3.48) 
Given the one-variable functions in (3.36), we have the following solution of (3.48): 

= 51 ^^^A'^)i€M'i-^'t)Ci,m,-i{a2,z)+^f3e,j{w)Qej^^^ (3.49) 

e=l,2 i=0 j=0 

In summary, we have the following theorem: 

Theorem 3.2. Let A he any function of t,z and let c be arbitrary real constant. Suppose 
the one-variable functions in (3.36) are given. We have the following solution of the nonlinear 
MHD equations (1.2) and (1.3): 

ip = atxy + Xtie'^x + ce'^y) + Ft{t, z, e^x + ce'"), (3.50) 

■0 = (Tzxy + Xzie^x + ce'^y) + Fz{t, z, e^x + ce'"), (3.51) 

where (1) a is given in (3.28) or (3.29) with ai,a2 G M and F{t,z,vu) is given in (3.37) via 
(3.32)-(3.34); (2) a is given in (3.38) or (3.39) with ai,a2 G M such that 02 ^ {0,Z + 1/2} 
and F(t,z,zu) is given in (3.44) via (3.32)-(3.34) and (3.41)-(3.42); (3) a is given in (3.46) or 
(3.41) with 01,02 G M\ G {0,Z + 1/2} and F{t,z,w) is given in (3.49) via (3.41) and (3.42). 

Finally, we have the following obvious results. 

Proposition 3.3. Let F{w, w) and G{w, w) be any two-variable functions. We have the 
following solutions of the nonlinear MHD equations (1.2) and (1-3): 

(1) if = Fy,{t -\- z,x) + Gw{t - z,x), ip = Fyj{t + z,x) - Gyj{t - z,x); (3.52) 

(2) 

(f = F^{t + z,x^ + y^) + G^t -z,x^ + y^), (3.53) 

^ = F^{t-{-z,x'^ + y2) -Gy,{t- z, + (3.54) 

We remark that we would get more sophisticated solutions if we apply the transformations 
(1.6)-(1.16) to our above solutions. For instance, applying the transformations (1.11)-(1.16) to 
the solution in (3.52), we obtain the following solution of the nonlinear MHD equations (1.2) 
and (1.3): 

(p = Fu,{t + z,ix + a) cos2a + {y + T)s'm2a) + a' i{x + af + {y + Tf ) 

+at{y + t) - Ttx + Xt + Gyjit - z, (x + a) cos 2a + (y + r) sin 2a), (3.55) 

tp = Faj{t + z,{x + a) cos2a-\- {y + T)sin2a) + ea'dx + af + {y + rf) 

+at{y + t) - Ttx + Xt - Gw{t - z, (x + a) cos 2a + (y + r) sin 2a) (3.56) 

with e = ±1, where a = a{t + ez) is an arbitrary one-variable function, and a,T,X are any 
functions in t,z. 
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